BONN-TH-2003-01 
FZJ-IKP-TH-2003-3 



Supernovae and Light Neutralinos: SN1987A Bounds on 
Supersymmetry Revisited 

H. K. Dreiner", C. Hanhart'', U. Langenfeld", and D. R. Phillips'^ 

"Physikalischcs Institut, Univcrsitat Bonn, NuBallec 12, 53113 Bonn, Germany; 
''Institut fiir Kcrnphysik, Forschungszentrum Jiilich, 52428 Jiilich, Germany; 
'^Department of Physics and Astronomy, Ohio University, Athens, OH 45701, USA; 



Abstract 

For non- universal gaugino masses, coUider experiments do not provide any lower bound 
on the mass of the lightest ncutralino. We review the supersymmetric parameter space 
which leads to light neutralinos, < 0(1 GeV), and find that such neutralinos arc almost 
pure bino. In light of this, we examine the neutralino lower mass bound obtained from 
supernova 1987A (SN1987A). We consider the production of binos in both electron-positron 
annihilation and nucleon-nucleon binostrahlung. For electron-positron annihilation, we take 
into account the radial and temporal dependence of the temperature and degeneracy of the 
supernova core. We also separately consider the Raffelt criterion and show that the two lead 
to consistent results. For the case of bino production in NN collisions, we use the Raffelt 
criterion and incorporate recent advances in the understanding of the strong-interaction part 
of the calculation in order to estimate the impact of bino radiation on the SN1987A neutrino 
signal. Considering these two bino production channels allows us to determine separate and 
combined limits on the neutralino mass as a function of the selectron and squark masses. For 
~ 100 MeV values of the selectron mass between 300 and 900 GeV are inconsistent with 
the supernova neutrino signal. On the other hand, in contrast to previous works, we find that 
SN1987A provides almost no bound on the squark masses: only a small window of values 
around 300 GeV can be excluded, and even then this window closes once > 20 MeV. 

1 Introduction 



If the supersymmetry-breaking, electroweak, gaugino masses fulfill the grand-unified mass rela- 
tion: 

Mi = -tanXM2, (1) 

then, because LEP chargino and neutralino-pair-production searches set a lower bound on M2, 
there is a concomitant lower bound on Mi. These constraints on the mass matrix result (for the 
case of conserved R-parity) in a limit on the mass of the lightest neutralino, x 

M^>46GeV. (2) 

However, Eq. is not an inescapable consequence of unification. For example, unification 
might occur through a string theory without a simple gauge group . Refs. 1^ showed that 
if the assumption ^ is dropped and Mi and M2 are considered as independent free parame- 
ters then there is currently no lower experimental bound on the neutralino mass from collider 
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experiments . In light of this situation, here we reconsider the bounds on neutrahno properties 
which can be determined from supernova 1987A (SN1987A) |Slini[2n8|. 

The basic idea is that in a supernova, neutrahnos with masses of order the supernova core 
temperature, = 0(30 MeV), can be produced in large numbers via electron-positron annihi- 
lation ISIIZI and nucleon-nucleon {NN) "neutralinostrahlung" jH]: 

e+ + e" — > X + X, (3) 
N + N ^ N + N + x + X- (4) 

Once produced, the neutrahnos have a mean- free-path, A;j^, in the supernova core which is 
determined via the cross sections for the processes El 13 ■ 

X + e — > x + e, (5) 
X + N ^ X + N, (6) 

as well as the electron and nucleon densities. If A-j^ is of order the core size, Rc = 0(10 km), 
or larger, the neutrahnos escape freely and thus cool the supernova rapidly. As the supernova 
temperature drops, the neutrino scattering cross section also drops and the neutrinos are no 
longer trapped, thereby further hastening the cooling of the supernova. Thus neutrahno cooling 
could significantly shorten the supernova neutrino signal 9 , in disagreement with the obser- 
vation from SN1987A by the Kamiokande and 1MB collaborations |10| lllj . However, if the 
neutrahnos have masses much greater than the supernova core temperature Tc then their 
production is Boltzmann-factor suppressed and they affect the cooling negligibly, independent 
of A;^. Demanding that be large enough that neutralino-cooling does not markedly alter the 
neutrino signal — particularly its time-structure — allows us to set a lower limit on the neutrahno 
mass. As we discuss in detail below, this limit depends on the squark and selectron masses 
which determine the relevant neutrahno cross sections for the production processes Q and (jlj, 
as well as those for the scattering processes © and 

A full treatment of this physics requires the implementation of the neutrahnos into the 
supernova simulation code. This is well beyond the scope of this paper. Instead we invoke 
two different analytic criteria to estimate whether neutrahno emission will affect the detected 
neutrino signal significantly. The first test we consider is to require the integrated supernova 
energy emitted through the neutrahno channel to be less than 

i.e. much smaller than the energy emitted by all neutrino species: Ei, ~ 3.0 x 10^^ erg The 
exact number we choose here is somewhat arbitrary-not least since it is subject to our ignorance 
of the total energy released in SN1987A. We will discuss the dependence of our neutrahno mass 
bounds on E^^^ below. 

A 

Detailed supernova simulations with axions suggest that another way to place bounds on 
these energy- loss mechanisms is the "Raffelt criterion" |12j . This states that "exotic" cooling 
processes, such as those considered here, will not alter the neutrino signal observably, provided 
that their emissivity, £, obeys: 

£ < 10^^ ergs/g/s. (8) 

If their emissivity is larger than this they will remove sufficient energy from the explosion 
to invalidate the current understanding of a Type-II supernova neutrino signal. With this 

^The sensitivity at LEP2 to the cross section a{ee — > X1X2) in the case where M^o and M^o are free parameters 
^ can be compared with the required sensitivity shown in Fig. 6 of Ref. 0. It is then clear that the light neutralino 
is not excluded by LEP2 but could possibly be produced at a next-generation collider. 
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simple criterion the results of the detailed simulations in Refs. |13[ll4j — where the exotic cooling 
mechanisms considered were axions |13| and KK-graviton emission |14| — can be reproduced. 
These detailed simulations also show that the feature of the SN1987A neutrino pulse which is 
most affected by the presence of alternative cooling mechanisms is the temporal distribution. 
Losing 10^^ ergs to neutralinos during the first ~ 1 second after bounce, as in our first criterion 
((7|). might not greatly affect neutrino temporal and spectral features. But, at later times, the 
loss of just a fraction of this energy could have a dramatic impact — by ending the period of 
diffusive neutrino cooling much earlier than would otherwise be the case. In this paper, we 
examine the impact of the process (jSj) using first Eq. (jT)) and then the Raffelt criterion, Eq. ©, 
and compare the results. We find good agreement between the two. 

Neither of these criteria are particularly suitable in the case that neutralinos have a small 
mean-free-path, A,^ < 0{Rc), since the neutralinos' contribution to the proto-neutron star cool- 
ing is then diffusive, rather than radiative. The treatment of this "trapped" regime is beyond the 
scope of this study. In Sec. El we make an estimate of the conditions for trapping and then refer 
the reader to the literature for bounds on SUSY parameters under these conditions d . In 
this work, we focus on estimating the emitted energy in neutralinos, assuming they are weakly 
coupled to matter, and then applying the criteria ((T)) and ((HJ. 

Previous authors EI considered the case of a stable neutralino and were thus restricted 
to the mass regions |15| ^: 

< 100 eV or M^> 500MeV, (9) 

in order to avoid an over-closed universe. The heavier neutralinos, > 500 MeV, are irrelevant 
for supernova and thus Refs. [SliniIZ| focused on a massless neutralino. In this work we consider 
the range: 

0<M^< 200 MeV . (10) 

We avoid over-closure of the universe by allowing for the possibility of a small amount of R- 
parity violation j20ll21j . and thus assume that the neutralinos, although stable on the time scale 
of the supernova and collider experiments, are not stable on cosmological timescales. For such 
neutralinos, which we call "quasi-stable", the mass restriction @ does not apply ^. 

To summarize: in this work we consider a general lower mass bound on quasi-stable neu- 
tralinos as a function of the squark and selectron masses. We assume the bounds from the relic 
density are avoided through small R-parity violating couplings. This is similar in outlook to 
Ref. [S] where, in light of the Karmen time anomaly j24j, the impact on SN1987A of a neu- 
tralino with the specific mass {M^ = 34 MeV) was examined. We extend this work to a general 
neutralino. We also differ from previous work in the following more technical points 

1. Throughout we consider a bino neutralino instead of a photino, since in Refs. O |3| it 
was shown that a light neutralino, < 5 GeV, must be dominantly bino in order to be 
consistent with LEP results. This leads to a substantially weaker effective coupling to 
nucleons and thus to weaker bounds on the supersymmetric-particle masses. 

■^The lower bound has recently been reconsidered |l(illl7ini^ in the light of new data from LEP and [g — 2)^. 
Associating the neutralino with the dark matter of our universe these works obtain > 0(5 GeV). The upper 
bound in Eq. Q was not considered in Refs. 1171 IT^ . although in fact sufficiently small neutralino masses 
are always allowed by the closure constraint, independent of the sfermion masses. However this light a neutralino 
would constitute a sizable hot dark matter component, whose presence structure-formation arguments strongly 
disfavour |19j . 

detailed analysis of cosmologically allowed neutralino lifetimes |21ll22ll^ will be given elsewhere. 



3 



2. For the case of electron-positron annihilation and the criterion Eq Q, we consider the 
radial and temporal dependence of the temperature and the degeneracy in the core j25j . 
This gives a more realistic estimate of this contribution. It also gives a stricter bound 
than some previous work [S] since the temperature in the outermost region of the star 
(r ~ ^Rc), where most of the neutralinos coming from annihilation are produced, is 
somewhat higher than the temperature considered by, for instance, Ref. However, it 
gives a weaker bound than in the earlier work by Ellis et al, Ref. [HI, where a significantly 
higher core temperature was chosen. Since all emissivities vary rapidly with temperature 
this results in a marked difference in the bounds on the supersymmetric-particle masses. 

3. The calculation of the A^A^ NNxx emissivity by Ellis et al. was based on the compu- 
tation of the rate for the related neutrinostrahlung process: A^A^ — > NNi^v jSj- For this 
Ellis et al. used the pioneering calculation of Priman and Maxwell, Ref. |26j . A recent 
model-independent treatment of NN dynamics in the production of axial radiation — e.g. 
by reactions such as A^A^ — > NNvD or A^A^ NNxx — suggests that Ref. overesti- 
mates supernova emissivities by about a factor of four [JJIEHlEn]- This loosens the bound 
previously obtained from SN1987A on the relevant supersymmetric-particle masses. 

4. Since previous authors only considered very light neutralinos, was neglected in com- 
puting the total neutralino emissivity. We calculate the full M^-dependence of In 
particular, for > Tc the mass dependence is very strong. Indeed, ultimately it is 
exactly this strong mass dependence which we use to derive a bound on My.. 

5. We also combine the emissivities from both electron-positron annihilation and nucleon- 
nucleon neutralinostrahlung to get a bound on supersymmetric-particle masses from both 
sources of neutralino radiation. 

The outline of the paper is as follows. In Sec. 121 we discuss the possibility of a light neutralino 
in the MSSM. In Sec. 01 we determine the bounds which can be found from electron-positron 
annihilation. This allows us to exclude certain regions of the joint parameter space of selectron 
and neutralino masses. In Sec. jll we determine the effective neutralino-nucleon coupling. In 
Sec. we use this to compute the general bound obtained from A^A^ NNxx ^ & function of 

and Mg. This analysis can be applied to the production of any particle pair which is axially 
coupled to nucleons. In Sec. (HI we discuss the effect that neutralino trapping due to matter 
and/or gravitational interactions has on these bounds. Finally, in Sec. [3 we combine the results 
for e~^e~ XX ™d A^A^ NNxx order to get overall information on the regions of MSSM 
parameter-space which are excluded, and offer our conclusions. 

2 A Light Neutralino in the MSSM 

Before we discuss in detail the bounds we obtain from SN1987A on a light neutralino, we 
consider whether such a neutralino is consistent within the MSSM and with existing collider 
data. In Ref. it was shown that there are no laboratory bounds on a neutralino with mass 
= 34MeV, provided it is dominantly bino — mainly because a bino neutralino does not 
couple to the Z'^-boson at tree-level. For a neutralino with mass < C(200MeV) — which 
we consider here — none of the bounds in Refs. |31 01 depend sensitively on the mass. We thus 
expect such a neutralino to be consistent with all terrestrial data ^. 

We have gone beyond Refs. |H1I1], and looked in more detail at where in the MSSM parameter 
space a neutralino with Mj^ < 200 MeV can be obtained. We have found new parameter regions 

*A detailed analysis can be found in Ref. |3U|. 
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Figure 1: The neutralino mass as a function of M2 for various values of Mi and with tan/3 = 10, 
/i = 300 GeV. The downward spikes correspond to true zeros of the neutralino mass. The bold 
horizontal line indicates the value = 200 MeV, below which we consider in this paper. To 



the right of the bold vertical line at about M2 
mass bound from LEP2. 



120 GeV the chargino mass satisfies the lower 



which are shown in Fig. ^ Note particularly that for a fixed value of Mi there can be more 
than one value of M2 leading to a specific M^. It is also clear that, as was pointed out in 
Ref. [SJ, neutralino zero modes are possible. These are the downward spikes in Fig. ^ There 
are in fact extensive regions in the Mi — M2 parameter space where a light neutralino can be 
obtained. However, if we take into account the lower mass bound on the chargino from LEP2 
(M2 > 120 GeV ), we see that for^ tan/3 = 10 and /i = 300 GeV a light neutralino is only 
obtained for the specific range 1 GeV < Mi < 1.5 GeV. (These parameter values will be shifted 
by the one- loop radiative corrections to the gaugino masses.) This involves a relative fine-tuning 
of 10~^ between Mi and M2. While this is perhaps aesthetically displeasing, it is by no means 
forbidden. 

In this parameter range, the neutralino is more than 98% bino. In the following we shall 
thus work with a pure bino. Our only free super symmetric parameters are then Mj^ and the 
selectron and squark masses. 



3 Electron-Positron Annihilation to Neutralinos 

We first discuss the case of free-streaming neutralinos. We defer the issue of trapped neutrali- 
nos to Sec. ini where we consider their scattering off electrons and nucleons in the supernova 
environment, as well as gravitational trapping. 

^tan/3 = Vi/v2 is the ratio of the vaccum expectation values of the two neutral CP-even Higgs bosons in the 
minimal supersymmetric standard model. ^ is the mixing parameter of the Higgs fields in the superpotential. 
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3.1 Emissivity for Free Streaming Neutralinos 

The neutralinos are produced via electron-positron annihilation 

e+{pi) + e-{p2) ^X{ki) + x{k2). (11) 

where we have indicated the four-momenta of the particles. The dominant process proceeds via 
t- and -u-channel selectron exchange. The energy thus emitted by the supernova per unit time 
and unit volume is the emissivity 

^(M^,re,r/)^^ = |^^L^/i/2(i?i + i?2)|Av|<T(e+ + e--X + x). (12) 

Here Ei + E2 = + E4 is the combined energy of the emitted neutralinos. The Fermi-Dirac 
distributions are given by 

where ^ii is the chemical potential and Tc the temperature in the supernova. In the following, we 
shall write rj = n/Tc for the degeneracy of the electrons. We have neglected the Pauli blocking 
of the final state neutralinos: (1 — /3)(1 — /4). In Eq. ((T^ . a{e'^ + e~ — > X + x) is the (free-space) 
cross section for the process (jllj) and |Av| is the absolute value of the relative M0ller velocity 



VM^l = V(V1-V2)2-(V1 X V2)2 (1 - COS 6) . (14) 

(v, are the velocities of the incoming electron and positron and 9 is the angle between them.) 
In the case where the selectron masses are degenerate, Mg^ = Mg^ = Mg the cross section is 
given by 

aie^+e ^x + x) = ^ j U -\ , (15) 

^ ^ 24cos4 ^^M| y s J ' ^ ^ 

where we have used rrig <^ s <^ M?. Here a is the fine structure constant, 9^, is the electroweak 
mixing angle and s is the center-of-mass energy squared. Replacing the bino-coupling by the 
photino coupling this agrees with Ref. '32[. 

The emissivity in Eq. (|12|) is a function of the neutralino mass via the cross section (|15|1. It 
is a function of the supernova temperature and rj via the Fermi-Dirac distributions (jl^j) . 

3.2 Total Emitted Energy 

In order to determine the total energy emitted in the neutralino channel we must integrate the 
emissivity over the time, to, during which neutralinos are emitted and over the volume of the 
supernova core 

E^{M^) = f'dt [ d\ £{M^,T,{r, t), ry(r, t)) ^ dt P{M^, t) . (16) 
Jo J Jo 

When performing these integrations care must be taken since the temperature and the degener- 
acy depend on the radius and on the time. We use the temperature distribution given in Fig. 1 
of Ref. |25j . There the temperature is shown as a function of the enclosed baryon mass. In Fig. 5 
of Ref. [ini we see that for t > 250 ms the density is constant. This makes the conversion of d^r 
to dM straightforward. Here we adopt a core radius of Rc = 13 km and a mass of M^tss = 1.4 Mq 
and obtain 

/> « 3 X 10^^ g/cm^ (17) 
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Figure 2: Neutralino power P^{My.,t) for t = 0.5 seconds as a function of the neutralino mass, 
My., with the selectron mass Mg = 200 GeV. The solid curve shows our numerical computation 
while the dashed curve represents the fit Eq. Hf^) . The two curves are almost indistinguishable. 



(Throughout we assume radial symmetry j25j.'l 

The time t = of Ref. ^5] corresponds to the time when the incoming shock wave stops 
and bounces outwards again. The boundary conditions at this time are not well known (33j . 
However, the exact shape of the initial conditions has little effect on the subsequent evolution 
and after 0.5 s the distributions are reliable. For the first second we have used the t = 0.5 s 
distributions of Ref. [221 • ^'^^ longer time scales we have used the subsequent radial distributions 
for Tc(r, t) and ry(r, i). However, as we will show below, most of the energy is emitted during 
the first second after the bounce and so ultimately we will use tQ = 1 second in (fTKI) to derive 
our neutralino mass bounds. This means that we are demanding that the prompt neutralino 
pulse in the first second does not have a total energy greater than the bound ®. We discuss 
below how the neutralino-mass bound depends on the value chosen for to. 



3.3 Results 

Before presenting the total emitted energy, we consider the power radiated in the neutralino 
channel. In Fig. 121 we show as the solid curve the power of the neutralino emission at t = 
0.5 s as a function of the neutralino mass, P(M^,0.5 s). We have fixed the selectron mass to 
Me = 200 GeV. According to Eq. ((T3|) . the power scales as Mg~^. A good fit to the solid curve 
throughout most of the neutralino mass range is given by 

„, , ^/200GeVV / 4 3 2 X My 

= ^ — ITJ exp(a4X + a^.x'' + a2X + aix) , x = , (18) 

\ Me / MeV 



ai = 1.71 X 10~^ 02 = -1.58 x 10"^ 



where the constants are given by 
1.71 X 10~^ , 

as = 2.695 x 10""^ , = -1.99 x 10"^° , C = 1.6318 x 10^^ erg/s. (19) 

We have compared our integration with the ratio presented in Ref. IBI and agree EH. We thank Michael 
Kachelriess for discussions on this point. 
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Figure 3: The dependence of the neutralino lower mass bound on the selectron mass. Here we 
chose E'^^^ = 10^^ erg and t^ = 1 s. The green-dashed line indicates the lower hound on the 
selectron mass from LEP2. Beyond Mg = 1 TeV the exclusion curve drops rapidly. 

The numerical fit in Eq. (|18j) is shown as the dashed curve in Fig. |2l Throughout the mass range 
considered it agrees to better than 1% and the two curves are almost indistinguishable. 

When using the Raffelt criterion it is also convenient to have a parameterization for the 
emissivity 

^(M^) = (^^^^^^) eMb4X^ + hx'' + b2X^ + bix), ^=^^ (20) 

and the constants are given by 

bi = 4.75 X 10"^ , b2 = -2.25 x 10"^ 

63 = 4.66 X 10"'^ , 64 = -4.02 X 10"^° , D = 9.0125 x 10^^ erg/g/s. (21) 

Once we have computed Ey.{My.) we can determine a smallest permitted neutralino mass by 
requiring the emitted energy to be below i??^^ of Eq. ((T)). If we choose i??^^ = lO^'^ erg, to = 1 s 
and the neutralino to be pure bino then for Mg = 300 GeV we have: 

^mm^ 243MeV, (22) 

while for Mg = 1 TeV we find M?™ = 90 MeV. For a massless neutralino to be allowed we 
require Mg > 1275 GeV. 

The value of M?™ is shown for a wide range of selectron masses in Fig. |31 Note that values 
of My. below the solid (red) line are forbidden, as are values of Mg to the left of the solid line, 
since in either case the total energy produced by the process in the first one second will be 
larger than i??^^. The lower selectron mass bound from LEP2 |35j : 

Me > 99.6 GeV, (23) 

is also indicated by the vertical dashed (green) line in the figure. 

How sensitive are these results to the choice to = 1 second? We made this choice because we 
expect most of the neutralino power to be emitted in a burst during the early, hottest, part of 



8 



300 



^00 



Si 00 



— M- 



300 GeV 



M-= 500 GeV 
-. M- = 1000 GeV 

e 



1,5 



2,5 



t (j(sec) 
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±y±^ , as a function of the time over which 
we integrate, to- The curves correspond to three different values of the selectron mass: Mg = 
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10 erg and the neutralino is pure bino. Our choice of to = 1 



second is indicated by the grey vertical line. 



the supernova. This is indeed the case. If we integrate the power out further, thereby increasing 
to, and again apply the criterion Q, we obtain the bounds shown in Fig. ^ Once tQ > 1 second 
the bound is essentially independent of to- 

As discussed above, the exact value chosen for E^'^^ — 10^^ ergs — is somewhat arbitrary. In 
Fig. IS] we show the dependence of the lowest- allowed My. on the choice of E^^^. This is done 
for several values of the selectron mass. Note that if E^^^ is decreased (increased) by a factor 
of two then M™° becomes at most 25% larger (40% smaller). 

Thus, the main dependence of the neutralino mass bound is on the selectron mass, which 
sets the effective electron-neutralino coupling strength. In addition, it is important to note that 
most of the neutralino production from electron-positron annihilation occurs in the outermost 
10% of the star {jqRc ^ ^ ^ Rc)- Looking at the profiles used for density and temperature 
during the first one second j25j the reason for this becomes apparent. At early times, this outer 
region has the highest temperature, and, more importantly, the lowest electron degeneracy, 77. 
Since the rate of the process Hll|) is suppressed by factors of e~'^, the smaller values of rj near the 
surface of the star mean that the majority of neutralinos are produced there. This will prove 
critical when we look at the effect of neutralino trapping on our bounds in Sec. 

3.4 Raffelt Criterion for e"'"e~ XX 

Although we have used the best supernova input available for radial temperature and density 
profiles, and have argued that our results are not especially sensitive to our choices for E^^^ and 
to, we would like to check the bounds on as a function of Mg that resulted from our modeling 
of the supernova. To this end we now turn to the Raffelt criterion 0. This is an estimate of 
what would happen were the neutralinos implemented in a full simulation of SN1987A. It is a 
test that is local in space and time — unlike the integral measure 0- 

The one parameter we are free to choose in applying Q is the temperature at which is to 
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Figure 5: The neutralino mass bound as a function of E^^^ for three different values of the 
selectron mass: Mg = 300, 500, 1000 GeV. Here to = 1 s and the neutralino is again pure bino. 
Our choice of Eq. ^ is represented by the black vertical line. 
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Figure 6: The neutralino lower mass bound Ss a function of the selectron mass Mg, when the 
Raffelt criterion is applied to the emissivity for electron-positron annihilation to neutralinos. 
The vertical line indicates the lower bound on the selectron mass from LEP2 1,^.5^ . 



be computed. Previous work suggests Tc = 30 MeV is a reasonable choice ^lEHSEl- Computing 
the emissivity due to e'^e~ — > xx this temperature, and at a density of 3 x 10^^ g/cm^, 
produces the constraint on neutrahno and selectron masses shown in Fig. IHl The solid (red) 
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Figure 7: Two diagrams contributing to neutralinostrahlung from the nucleon-nucleon system. 
Time flows from left to right. The nucleons are denoted by double lines. The (yellow) oval 
indicates the nucleon-nucleon interaction. The black circle indicates the effective neutralino- 
nucleon coupling which we compute in Sec. ^ Two diagrams which are the same, except that 
the neutralinos are radiated from the bottom lines, are not shown. 

line indicates the neutralino and selectron masses for which the emissivity of the process is 
exactly 10^^ ergs/g/s. Neutralino masses below the solid (red) line, and selectron masses to the 
left of it, are forbidden. According to the criterion ((HJ, if the supersymmetric particles had such 
masses the time-structure of the observed SN1987A neutrino pulse would have been noticeably 
different. 

Note that the bound obtained with the Raffelt criterion is a little less stringent than that 
obtained from ^ and depicted in Fig. |3I The two criteria would be in agreement if we chose 
^max roughly a factor of two larger. Alternatively, if in the Raffelt criterion we chose the 

supernova core temperature to be T^ = 34MeV, we would also obtain agreement between the 
two approaches. 

Thus the qualitative agreement between these "global" and "local" criteria gives us confi- 
dence in the computationally-simpler Raffelt criterion — confidence that is supported by detailed 
simulations of the impact of axions JH] and KK-gravitons on the SN1987A neutrino signal. 
From now on we will use the Raffelt criterion to assess the impact of neutralinos on the cooling 
of SN1987A. 

4 Neutralino-Nucleon Interactions 

We now turn our attention to the "binostrahlung" process ^ . We wish to evaluate the emissivity 
due to this reaction. To do this, in SecElwe compute the Feynman graphs shown in Fig. [3 This 
gives the amplitude for the process (jlj) to leading-order in the soft-radiation approximation. 

These two diagrams, together with their partners under 1 <-> 2 interchange, give the dominant 
result for binostrahlung in the limit that the total neutralino energy a; ^ 0, since they are the 
only graphs which diverge 1/"^) in that limit. Indeed, were the radiation a single photon, 
rather than a XX pair, the graphs in Fig. [3 would yield the leading-order Low energy theorem 
for NN NN'j |27l 1371 IBSj . In order to evaluate them we need to know two amplitudes: the 
nucleon-nucleon interaction, denoted by a (yellow) oval in Fig. [3 and the effective neutralino- 
nucleon interaction, represented by a large black dot in Fig. [Tj 

In the soft-radiation limit the intermediate-state nucleon line in Fig. [3 is almost on-shell. 
Thus, as we explain in Sec. [3 and Appendix A below, to leading order in the soft-radiation 
expansion, the nucleon-nucleon interaction occurs on-mass-shell. In consequence, it can be 
reconstructed from A^A^ scattering data, giving a model-independent result for this piece of 
input to the diagrams in Fig. [3 
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Figure 8: The two diagrams contributing to neutralinostrahlung from a nucleon. The (green) 
oval denotes the nucleon bound state. The horizontal lines are the partons with q = u, d, s. The 
virtual squark q can be either left- or right-handed. 

In this section, we focus on the other piece of input to these diagrams: the neutrahno-nucleon 
couphng. We will show that in the limit of long-wavelength (i.e. soft) neutralinos this coupling 
is expressible in terms of the spin content of the nucleon, Aq. We evaluate the xn coupling using 
two different sets of results for An, Ad, and As: the predictions of the non-relativistic quark 
model (NRQM), and those of the quark-parton picture at leading order in (LO-QPM). The 
neutralino- neutron coupling turns out to be very different in these two cases, and this difference 
has a considerable impact on our final predictions for the emissivity due to nn binostrahlung. 
We believe this gives an estimarte of the range of permitted values. 

4.1 The Bino-Quark Four-Fermion Effective Lagrangian 

The Lagrangian for a pure bino-neutralino coupling to a quark and a squark is given by j39j : 
A^-x = -^^^—Y,[YgAQRPLX)qR + Y,,qmPLq)], (24) 



where e denotes the electric charge. The sum runs in our case over the quarks q = u,d,s 
and (Jl^r denotes the corresponding left- or right-handed squark. The hypercharges Yg^ for the 
right-handed singlet quarks are given by: ~ ~|' Ydj^ = Ys^ = |. For the left-handed 
quarks Yg^ = |. For each term in the Lagrangian we obtain two crossed Feynman diagrams 
contributing to neutralino pair emission from a nucleon, as shown in Fig. |H1 In the following we 
shall assume that the squark masses are degenerate: 

Mg = M^^ = M^^ = Ms, = Mu^ . (25) 

Since the momentum transfer in the neutralinostrahlung in our case satisfies <C M|, we can 
obtain an effective four-fermion operator for the qqxx interaction for a pure bino |4nj : 

= 8 cofe^ M? ^^^'^'^^ ^ ^^"^^^^^^ ~ ^^"^^^^^^ ] • ^^^^ 

Here we have performed a Fierz transformation in order to separate the neutralino current from 
the quark current. The vector current for the neutralinos vanishes since they are Major ana 
fermions. The quark current can be split into a vector and an axial-vector part: 

9 a 
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with the effective charges {W^ = _ y2 ^ = y2 ^ y2 



^u=f; W!=W^ = l, (28) 



^^=f; Xj=Xf = ^. (29) 

4.2 Nucleon Matrix Elements 

The neutrahno coupling to nucleons can be calculated by taking matrix elements of the inter- 
action in Eq. H27|) between nucleon states. In the limit that <^ ^qcd ^^^^ interaction energy 
can be rewritten in terms of an axial neutralino current coupling to quark vector and axial- vector 
current matrix elements: 

9 q 

with: 



Vi'^ = {n\ J d^xq{x)^^q{x)\n), (31) 
4") = {n\ [ d'xq{x)j^j5qix)\n). (32) 



These quark-current matrix elements simplify to: 

^;l^) = H^|n)7Vf), (33) 
4) = (n|Sjn)A(z("), (34) 

(n) 

where and M„ are the momentum and the mass of the neutron, A'^^ is the number of valence 
quarks of type q in the neutron, is twice the Pauli-Lubanski spin vector, and Ag^") is the 
average contribution of quarks of type q to the total neutron spin. 

For a neutron at rest, the only component of V^"^^ which is non-zero is the zeroth one, while 
the four-vector has only non-zero spatial components. Thus we arrive at an interaction 
energy: 



V< 



e2 



16M? cos2 Ow 



X7°75X ^ - ^775X • (nV„n) ^ Aq^^^ 



(35) 

with (Jn the three- vector of the neutron Pauli spin matrices. Corrections to Eq. pSj) are sup- 
pressed by powers of P„/-M„ and/or powers of QRn, with Rn the neutron's size. This allows us to 
write an effective four-fermion Lagrangian for the interaction of neutralinos with non-relativistic 
neutrons: 

(c^'5/x,o - c^Sfi,iCri) n , (36) 



r-- - ^ SUSY a 



with the axial- vector neutralino current: 

a^^lri^X, (37) 
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and the effective couplings: 



e2 



^v = f; 4 = ^Ad+^iAu+ As). (39) 

(In (|39p we have used the conventional definitions of the Ag's and isospin.) 

The two diagrams in Fig.EJenter with opposite signs in the computation of the amplitude for 
the radiative process. Thus, as we discuss further in Sec. El to leading order in 1/Mn the vector- 
current interaction does not contribute to the neutralinostrahlung emissivity. Hereafter we 
disregard its contribution in order to get the lowest-order result for the binostrahlung emissivity 
in the soft-radiation approximation. It is worth bearing in mind though that 1/M„ corrections, 
as well as effects due to non-zero proton fraction, could render vector-current radiation quite 
significant in binostrahlung. 

Note that, in general, the vector current dominates neutralino- neutron scattering, which we 
will examine in Sec. El 

This leaves us with the task of estimating Au, Ad, and As. We consider only two out of a 
plethora of approaches: 



1. Employ the non-relativistic quark model (NRQM). In this case: 

-; Ad = — ; 
3 3 



4 1 

Au = -; Ad = — ; As = 0, (40) 



giving: 
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(41) 



2. Employ the leading order quark parton model (LO-QPM) j^. This obeys the constraint 
Au — Ad = QA, and does not fare too badly with respect to the constraint from the F/D 
ratio extracted from hyperon decays i42j: 



Am + Ad - 2As « 0.682 . (42) 

The values from Ref. |1J are: 

Ati = 0.78 ±0.03; A(i = -0.48 ± 0.03 ; As = -0.14 ± 0.03 . (43) 

These values rely on the extraction of the flavour-singlet matrix element from the zeroth 
moment of the protons gi structure function, and its subsequent evolution to low-Q^ . This 
evolution is perhaps questionable, since below ~ 1 GeV^ it becomes non-perturbative. 
Nevertheless, the values give: 

= -0.55 ± 0.06 , (44) 
a result five times as large as the NRQM one, and of opposite sign. 

4.3 Ratio of Neutralino and Neutrino Processes 

We now follow Ref. 0, and take the ratio of the coupling Ca to the corresponding neutrino 
coupling, in an attempt to eliminate some of the uncertainties in the computation. 
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Ca 



^ NRQM, 



The effective Lagrangian for tlie neutrino case is 1271 143j : 

G p -I- 

with the leptonic current: 

= 9-i^'{l--i^)v . (46) 

Here we will drop the vector-current piece, since it does not contribute to nn nnvu. Mean- 
while, arguments analogous to those of the previous subsection yield: 

< = -[Au- Ad + As]. (47) 

In the non-relativistic quark model we have: 

< = -^, (48) 

while using the values for the Ag 's from the quark-parton model in Eq. (|43)1 we find: 

< = -1.12 ±0.05. (49) 

The latter value is about 33% lower than in the NRQM: a significantly smaller discrepancy than 
in the neutralino case. 

Therefore for our chosen evaluation procedures the ratio of neutralino to neutrino axial 
coupling constants is: 

rX 

(50) 

0.49 ±0.08 LO-QPM. 

In order to use the the effective Lagrangian H45|) to compute the neutrino emissivity due to 
the reaction nn — > nnvv, we require the spin-summed matrix element for nn nnvv |26 | I27 | H3|: 

\M? = ^^Tm,){^)m,. (51) 

u spin 

Here Tiij is the hadronic tensor, which corresponds to two insertions of the neutron spin in the 
hadronic (in this case nn) state. The evaluation of this quantity will be discussed in Sec. El but 
note that we have already used the fact that only its spatial components are non-zero in the 
non-relativistic limit. (Here and below we denote spatial indices by Roman letters.) Also, for 
later convenience, we have removed a factor of l/o;^, to take care of the leading behavior of |A^p 
in the soft-radiation limit, as well as a factor of 16, to account for identical particles and the 
fact that the neutrinos couple to both nucleons. {uj is the total emitted neutralino energy.) The 
sum over hadronic states in Tiij will then only include states permitted by the Pauli principle. 
Meanwhile, TT(lilj) are the space-space components of the neutrino tensor: 

i''" = Tr[/^i7^(l - 75) fe7.(l - 75)] , (52) 
= 8 [fc'/A:^ ± k!^k'( - g'^'ki ■ k2 + ie^^p^kfk^] , (53) 

in agreement with Friman and Maxwell [5^ (^1,2 are the outgoing neutrino momenta). 

Now, in order to compute the spin-squared matrix element for nn nnxx we proceed by 
analogy with this neutrino calculation. The spin-summed matrix element for neutralinostrahlung 
which is integrated to yield the neutralino emissivity due to the process (UJ) can be written as: 

^ |_^|2 ^ ^k^T.{a,a,) . (54) 

X spin 
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Crucially, the hadronic tensor here is the same as the one in Eq. (|51|1 . since neutralino-pair 
production is also induced by two insertions of the neutron-spin operator in the hadronic system. 

We now compute the neutralino tensor Tr(ajaj). It is to be understood as the space-space 
piece of the four-tensor: 

M^' = Tr[(^i + M^)7^75(^2-M^)7^75], (55) 
= 4 [k'^kl + k\ki^ - {ki ■ k2 - Ml)g^'^] , (56) 

where ki and k2 are now the outgoing neutralino four-momenta. 

In order to compare the neutralino and neutrino computations it is convenient to define 
to be the ratio of neutralino to neutrino spin-summed matrix elements, with the ratio taken in 
the hmit My. = 0. Then: 

2 



where Tiij is the hadronic spin response in both cases. 

In the computation that follows in Sec. El we assume that the nucleons are in thermal 
equilibrium, with a common temperature Tc- Thus the nucleon kinetic energies are of order 
Tc, and their momenta are of order ^/MnTc- The centre-of-mass energy of the reaction Q 
is therefore also thermally distributed around 2Tc. This means that neutralinos are produced 
with a common "temperature" ~ Tc- For My- < 50MeV, the neutralinos are then relativistic 
and have momenta of order <C ^M^T^. Consequently, to leading order in \jTcjMn we can 
ignore the recoil of the hadronic system when the neutralinos are emitted, and Hp.y will then 
be independent of ki and k2. This allows us to perform the integral over these variables and 
obtain an "angle-averaged" neutralino tensor: 

d?Lx f dVL2 
4tt / 4tt 



AA,. ^ — ^ -^M,. , (58) 



and analogously for the neutrino case. For the space-space components — which are all that are 
relevant here — this gives: 



Afij = 6ijA{ujiu;2 - M. 



6ij8uJiLL!2 , (59) 



where uji^2 are the energies of the two emitted neutrinos/neutralinos. 

In any observable — including emissivities — in which the outgoing neutr(al)inos are not de- 
tected we can use these "angle-averaged" leptonic tensors in the evaluation of the spin-summed 
matrix element — provided, of course, that we truly are in the kinematic regime where recoil of 
the hadronic system can be neglected. Calculating a^- using the tensors (ISU)) we get: 

IfGsusYclY ^g^^ 

(61) 



2V G^c- 
sin^ 9w fMzY fc^ 



M^ 



Note that the hadronic tensor has canceled, so this ratio is independent of the model of the 
nuclear medium that is adopted. Note also that all of the nuc/eon-model dependence is now 
isolated in the last factor. Finally, Eq. (|6T|) gives: 



(M^ 



4 { ^ NRQM, 



sin'^H/l^l x<' (62) 
I i LO-QPM, 
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where we have written 0.49 ^. is more than fifty times larger in the LO-QPM than in the 
NRQM. 

5 Emission from a Neutron Gas 

We are now ready to compute the emissivity due to neutron-neutron "neutrahnostrahlung" , the 
process @. (For simphcity we wiU restrict ourselves to pure neutron matter. The generalization 
to arbitrary proton fraction is straightforward, if somewhat cumbersome.) We shall derive a 
general emissivity, which is valid for any particle-pair which couples axially to the nucleons. We 
then demand that this emissivity is not larger than the Raffelt bound (jH]), and so set bounds on 
a^. At the end of this section we use the results of Sec. I^to translate our bound on into a 
bound on the squark mass. 

It was recently shown that, for soft radiation, the neutrino emissivity due to nn — > nnvv 
can be directly related to the on-shell nucleon-nucleon data [2Z1- For < 50MeV, this "soft- 
radiation approximation" also provides a model-independent way to compute the emissivity due 
to neutrahnostrahlung. On the other hand, for neutralino masses of order 100 MeV or more the 
radiation cannot be regarded as soft, since ~ m-,^^ and the pion plays an important role in 
A^A^ scattering and in bremsstrahlung dynamics. However, even for large (from a nuclear physics 
point of view) neutralino masses we still regard the method presented here as a useful way to get 
an order of magnitude estimate of the emissivity due to nn — > nnxx- As will be demonstrated 
below, this emissivity depends sensitively on through the Boltzmann factors and so, even if 
our estimate of the matrix element is off by an order of magnitude, the neutralino-mass bound 
changes by only about 30%. 

The emissivity for the radiation of neutralinos from an interacting neutron gas is given by 

(27r)V4)(pi+P2-yi-p'2-fe) 

X S ■ f{E,)f{E2){l - f{E[)){l - f{E'^))^ , (63) 

where S = 1/8 is the symmetry factor taking into account that the initial and final nucleon 
pairs, as well as the neutralinos are identical, pi, p2 are the incoming andp'^, p'2 are the outgoing 
neutron four-momenta. f{E) is the Fermi-Dirac distribution function (|13j) . The energy per unit 
frequency interval for a particular two-neutron phase-space element is given by: 

oj and k are the total energy and three-momentum of the neutralino pair, and ki^2 are the 
three- momenta of the neutralinos. M is the matrix element for the process nn — > nnxx- 

To proceed, we now use the decomposition of into a hadronic tensor and a neutralino 
tensor as given in Eq. H54|) . Introducing this into Eq. H64() and using the angle-averaged 
expression of Eq. ((S5|l for the neutralino tensor, we derive 

^ = ayGpc'^'^^ f duJiduJ2 5{uj - wi - ^2) (^1^2 - M?)7iii , (65) 

aU! J UJ ^ 

Here a^, Gp, c^^ are as defined in Sec. 01 We will now set bounds on which apply to any 
axially-coupled light-particle-pairs emitted from the neutron gas. We then translate to bounds 
on the squark mass Mq using Eq. (|62j) . 



£ 



doj 



n 



=1,2 



d^Pi S'p\ 
(27r)6 



17 



5.1 Evaluating the Hadronic Tensor 

In Section \^ we saw that when calculating the neutr(al)ino emissivity only the space-space 
diagonal elements of the hadronic tensor are required. To leading order in the soft-radiation 
approximation, and provided many-body effects are not large, only the diagrams shown in 
Fig. [7| contribute to Tin. The tensor can then be expressed in terms of the commutator of the 
production operator with the NN scattering matrix Tj^n [13 • According to Cj^xnn of Eq. 
in the non-relativistic limit only two operator structures govern the XX pair's coupling: the unit 
operator and the spin operator. The commutator of the unit operator with T/vat vanishes, and 
therefore the vector current in Eq. (|H6j) does not contribute to nn — > nnxx- Thus, at leading 
order the only structure which contributes to TCij is the two-nucleon spin operator, Si. For a 
given two-nucleon state of total energy E, we find that the relevant, diagonal, components of 
Hij are: 

nu{E)= \{lm'„p'\[S^,TNN]\p,lm,)W (66) 

where |p, Inis) is a two-nucleon state of total spin 1, total spin projection m^, and relative 
momentum p. Here, because the binos radiated in Fig. [7| are soft, the initial and final nn 
relative momenta obey: 

p2 = p'2 = ME . (67) 

The key point here is that the A^A^ matrix elements in Eq. (|66|) are on-shell, and so can be 
obtained from nucleon-nucleon scattering data. So, to leading order in the soft-radiation ap- 
proximation, there is a direct connection between this data and the neutralino emissivity d£/dt: 
a connection given by Eqs. and 

To further simplify Eq. (|63j) . spherical symmetry and energy- momentum conservation can be 
exploited, thereby eliminating nine of the eighteen integrals. We define total, initial relative, and 
final relative momenta P, p, and p', respectively, and also move to dimensionless variables |44j : 

T = , = , a = , (68 

8MTc 2MTc 2MTc ^ ' 

and 

r/n = ^, M^ = ^, (69) 

where Tc is the temperature of the neutron gas in the supernova. At large temperatures, the 
angular dependence of the Fermi functions becomes weak — in the non-degenerate limit the 
angular dependence disappears completely. We have checked that at the temperatures we are 
studying replacing the product of Fermi functions by their angular average is justified. These 
angular integrals can be performed analytically 



' :t: f(F )f(F)- In f + + ^^^^ 1 (70) 

where = i((5-)-(j)2_j_-j-2_^^_ Doing the angular integrals this way necessitates the introduction 
of the angle-averaged hadronic tensor: 

/dp dp^ 
16^ ■ 

The replacement Tin — > Tin is not only a good approximation, it is a useful one, since the 
integrals in Eq. ()71() be carried out analytically, when the A^A^ T-matrix is given in partial wave 
form. The procedure used to construct Tin from a given set of A^A^ phase shifts is sketched in 
Appendix A. The results that follow were generated using the VPI phase-shift analysis in that 
construction [IH] , 
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Figure 9: Neutralino emissivity in units of 10^^ ergs/g/s as a function of the neutralino mass 
M^. The red solid line shows the full result for a^- = 10 a'^^^^, whereas the black dashed line shows 
the fit according to Eq. ( |77| j. The horizontal blue solid line shows the emissivity 10^^ ergs/g/s 
that is used when the Raffelt criterion is applied. 



5.2 Emissivity 

We now return to the computation of the emissivity Eq. (|63|) . Doing the angular integrals as 
described leaves us with a four-dimensional integral which we decompose as 



8(2 



^h{5/M^ 



da F{5, fj, rjn ) Ha {Tc a) 



where 



with 



h{y) = -{y- If + l)I,(y) _ (y _ lfl,{y)] , 



h{y) 



and Yiy) = (y + l)/(y 
T{5, a, r]n 



1). The result of the angular integration of the Fermi functions is 



drTZ 



where 



n{a, T, ri„ 



-{a + 5),T,r], 



In 




5),T,ri„ 



( cosh {\{(T + af - r]n)) 



sinh(a2 + - I cosh (^((r 



Vn)) 



(72) 

(73) 
(74) 

(75) 
(76) 



The functions Ik{y) can be evaluated to very high accuracy if we replace the second factor in 
the square root by its Taylor expansion up to third order in {x/Yf. Formally one would expect 
this approximation to work for Y ^ 1 only. However, even for Y = 1 both I2 as well as Iq turn 
out to be within 5% of the exact result. The remaining 3 integrals are evaluated numerically 
using the Gauss-Legendre method. The computation was performed for nuclear matter density 
and Tc = 30 MeV. The result of the numerical integration is shown as the solid line in Fig. 121 for 
Oy = 2.5 X 10~^. The numerical result is well fit by: 
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Figure 10: The region of the a^-My.(MeV) plane that is excluded by our soft-radiation calcu- 
lation, when applying the Raffelt criterion (at = 30 MeV ) is indicated by the grey shading. 
The blue solid line is the curve that defines the values of and which produce a neutralino 
emissivity of exactly 10^^ ergs/g/s. The horizontal (red) dashed line indicates the uppermost 
value of at which the soft-radiation approximation should be trusted. 

£ = 4a^ {I + 2M^ + 0.55M?) e'^^^^ x lO^^ [ergs/g/s], (77) 

where = M^/T^ as defined in Eq. The fit is good for < 200 MeV. Our result 

represented in the form (|77() allows us to directly derive a neutralino mass bound for any given 
coupling strength using the Raffelt criterion ©. Alternatively, we can think of defining a 
function a"'*(Mj^) from Eq. 1)77(1 : the coupling at which the Raffelt criterion is exactly met for 
a given neutralino mass M^. 

By setting Af^ = in Eq. ()77() above, i.e. a massless neutralino, we obtain the minimal 
coupling a^"*(0), for which we can set a bound. Using Eq. Q we obtain 

af\0) = 2.5 X 10-^ (78) 

For M^ = Tc = 30 MeV, af'\M^) increases to 5.2 x 10"^. For > Tc, the minimum we 
are sensitive to changes rapidly with the neutralino mass due to the exponential dependence 
seen in Eq. (|77|). Using the Raffelt criterion we can exclude a region in the a^~My. plane which 
is to the right of the curve defined by the function a'^^^{Mj^); this region is shown in Fig. 1101 

As stated in the beginning of this section, we believe the soft-radiation approximation we 
make in our computation is under control for < 50 MeV. This limiting mass value is indi- 
cated by the horizontal (red) dashed line in Fig. 1101 At this mass value, the smallest value 
we are sensitive to is 

af*(50MeV) = 1.2 X 10"^. (79) 

For > a™* (50 MeV), we can still use the soft-radiation calculation to obtain a reliable bound, 
as long as < 50 MeV. This excludes the shaded area in Fig. 1101 

For larger neutralino masses, the soft-radiation approximation breaks down, since the emitted 
energy is of order 100 MeV. An approximate bound for > 50 MeV can still be derived, be- 
cause most of the dependence of the emissivity comes from the phase-space integrals in Eq. ()63(1 , 
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and so using the form (|54|) to fix tlie nn — > nnxx matrix element sliould give results accurate to 
an order of magnitude, even if is as large as 200 MeV. Thus we argue that the entire region 
of — My. parameter space to the right of the solid (blue) curve in Fig. is also excluded by 
the arguments of this section, although less rigorously so than the grey shaded area. 

5.3 Bound on Squark Masses 



Using Eq. (|62|) to relate a^^ to Mq we can turn the bounds on we have obtained into bounds 
on the squark mass. 

1/4 ( 132 GeV, NRQM, 



2.5 X 10-5 

I 360 GeV, LO-QPM, 



where all the dependence on is now in the value of oi'^^^^ which can be easily obtained from 
Eq. H77|) and the Raffelt criterion. 

Eq. (|80|1 is significantly weaker than the bounds found in Ref. [£]. There are several reasons 
for this: 

1. As noted in the Introduction, Ellis et al. use a supernova core temperature of Tc = 
70 MeV — more than a factor of two larger than the one employed here. Since the emissivity 
(|77|) varies as a high power of the temperature this is a crucial difference. 

2. The axial hadronic response to which the neutralinos couple is about four times larger 
in the Ellis et al. result than in our improved computation, because Ref. ^ (effectively) 
employed the single-pion-exchange approximation to compute this response. 

3. Ellis et al. assumed a neutralino which was pure photino. The effective coupling for a 
pure-photino neutralino compared to a pure-bino neutralino is a factor 



2 



256 cos^ 0^ 1 1 « 10.28 (81) 



^Ad+ |Ati+ |As 



larger. (Here we computed analogously to the manner used to find in Eq. and 
used the LO-QPM values for the Ag's.) 

Our bounds on the squark mass are thus significantly less restrictive than those we found 
for the selectron mass in Sec. |21 



6 Neutralino trapping 
6.1 Diffusive Trapping 

The arguments of this paper apply only if neutralinos free-stream out of the supernova. In order 
to determine whether this occurs we must estimate their mean-free path 

Ax = ^ , (82) 

via the scattering processes given in Eqs. ()5I6() . Here n is the target particle number density 
in the supernova. For the electrons, we take ng = 8.7 x 10^^ m"^ and for the neutrons nn = 
5.7 X 10"^^ m-3 [^. In both cases we approximate the density to be independent of radius. 
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We shall first discuss the scattering on electrons and then estimate the case for scattering 
on nucleons. In both cases we use the optical depth criterion |25j : 

to determine whether neutralinos produced at a depth ro, free-stream out of the supernova or 
not. 

Any neutralino propagating out of the proto-neutron star can undergo both electron and 
neutron scattering. There are, however, two distinct classes of neutralinos, each of which is 
(mainly) produced in a different proto-neutron star region. The two classes also have different 
energy distributions, which is important, since the mean-free path depends via the cross section 
on the neutralino energy. Strictly speaking, we should average over the entire spectrum of pro- 
duced neutralino energies, but here we estimate the effect of scattering by fixing the neutralino 
energy at the peak of the energy distribution and evaluating cross sections, and hence mean- free 
paths, at that value. 

Thus for both electron and neutron scattering we apply the test (|83|) separately for two 
classes of neutralino: 

1. Neutralinos produced by e'^e~ XX- Sec. 01 we found that most of these are born in 
the outermost 10% of the supernova. Thus in the case of these "annihilation neutralinos" 
we take ro = 0.9 Rc in Eq. (jSHJ- As for their effective energy, in e^e~ — > xX; we must 
average over the thermal electron energies to obtain the neutralino energy. We estimate 
this as follows. The neutralino production cross section is proportional to Combining 
this with the energy dependence of phase space and the Boltzmann suppression 
exp{—E^/Tc) we expect a maximum number of neutralino pairs with a combined energy 
{E^^) = 4:Tc. Thus for the average individual annihilation neutralino energy we adopt 
(E^) « 2Tc. (For > 2Tc we just set = M^.) 

2. Neutralinos produced by binostrahlung. In this case the neutralinos are predominantly 
produced in the centre of the supernova and so we apply Eq. (|83j) with ro = 0. By a similar 
argument to that of the previous paragraph, we estimate {E^^^) = 5Tc, and thus for the 
individual neutralinos (E^) = ^Tc. 

In both cases the produced neutralinos can scatter off either electrons or nucleons. Thus, for 
each type of neutralino we now determine separate bounds from ex and nx scattering on Mg 
and Mq respectively. 

6.1.1 Neutralino-EIectron Scattering 

The scattering cross section for a pure bino on electrons is given by"^: 

''<' + '<^' + >" = iw«:m| ('-^j (' + ^j • 

where s = M| + 2EeiE^ - \p^\cos9), is the centre-of-mass energy squared. The neutralino 
mean-free path A-,^ depends strongly on the selectron mass, Mg, and, less strongly, on M^^. For 
a small selectron mass, we obtain a large cross section and thus a small mean-free path. In the 

^We obtain the scattering for a massless photino by setting = and multiplying by 2 cos* (ez,)* + 

y(efl)''), where Y{eL) = -1/2 and Y{eR) = -1. The resuH agrees with Ref. Q. 
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following, for a given neutralino mass, we shall determine a lower bound on the selectron mass 
above which free-streaming occurs. 

Estimating the electron energy as E^, = |T, we find that for free-streaming of annihilation 
neutralinos to occur, the selectron mass must obey: 



M^(MeV) 





10 


30 


50 


100 


150 
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Me(GeV) > 


322 


322 


321 


315 


341 


363 


377 



where we have written the results for values of between and 200 MeV. The fact that the 
bound first drops and then rises with rising M-^ is due to the functional dependence of the cross 
section ((SI)) . 

Me must be larger than about 320 — 380 GeV in order to ensure free-streaming. Thus 
even though Sec. |21 suggested that selectron masses of order 200 GeV were forbidden by the 
SN1987A data (provided that < 200 MeV) this is not strictly speaking the case, since if the 
selectron mass were that low then the neutralinos produced by electron-positron annihilation 
would propagate diffusively out of the supernova. To get an accurate bound for Mg ~ 200 GeV, 
neutralinos would have to be coupled to matter and included in a supernova simulation. 

There is thus a window of allowed selectron masses 100 GeV ^ Mg ^ 300 GeV where inter- 
actions with electrons mean that the annihilation neutralinos will diffuse out of the supernova. 
This physics is beyond the scope of this paper and we have nothing to add to the previous 
approximate computations of this regime performed in Refs. [SlISI- 

Meanwhile, a neutralino produced via binostrahlung can also scatter off the supernova elec- 
trons. These neutralinos free stream as long as: 



M^(MeV) 





10 


30 


50 


100 


150 


200 


Me (GeV) > 


606 


606 


604 


596 


607 


646 


670 



However, if Mg is lower than 600 GeV copious numbers of neutralinos will be produced 
by annihilation, and the diffusion of binostrahlung neutralinos can be safely neglected. So the 
possibility of an Mg low enough to induce binostrahlung neutralino diffusion is already excluded 
by our other arguments, unless M^ > 150 MeV. The bounds of Sec. 01 do not preclude Mg ~ 
650 GeV for this heavy a neutralino, but if > 150 MeV we will see below that binostrahlung 
cannot be used to set a bound on the squark mass anyway. Thus the constraint on Mg from the 
consideration of binostrahlung neutralinos is essentially irrelevant. 

6.1.2 Neutralino-Neutron Scattering 

The neutralinos can also be trapped via neutralino- nucleon scattering, Eq. @. In order to 
compute the cross section, we can employ the effective Lagrangian Eq. We obtain 

a{x + n^X + n) = [cf (S? - M?) + 3cf (^ + M|)] . (85) 

We have made the approximation for the centre-of-mass energy squared s « M^, where M„ 
is the neutron mass, and assumed non-relativistic neutron spinors. We see that we obtain 
contributions from both the vector and axial-vector neutralino current. The vector-current con- 
tribution vanishes in the limit of non-relativistic neutralinos — a limit that is taken, for example, 
in computations relevant to dark-matter detection. 
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For the case of neutralinos produced via neutralinostrahlung, we take {E^) = ^Tc. Setting 
ro = and using the results from the LO-QPM, we then obtain the lower limit on the squark 
mass 



M^(MeV) 





10 


30 


50 


100 


150 


200 


Mg-(GeV) > 


306 


305 


298 


280 


269 


330 


381 



Comparing this with the results in Sec. |SJ we see that even for the case of most dramatic 
neutralino-cooling effects — = — we can only exclude a region 300 GeV < Mq < 350 GeV by 
the arguments presented there. For > 30MeV, the values of Mg excluded by the supernova- 
cooling analysis of Sec. El are so low that the supernova neutralinos would be trapped by their 
strong interactions with neutrons in the proto- neutron star. 

It is interesting to note that scattering is competitive with production — even for = — 
partly because of the sizeable ratio of vector to axial-vector coupling constants. The large 
value of Cv/ca also means that if < fTc the cross section ()85() is dominated by the term 
proportional to Cv^. This piece drops with M^, but once > ^Tc the cross section begins to 
rise with M^^ as the axial contribution takes over. This interplay of vector and axial-vector pieces 
produces the dependence of the Mg-bound on seen in the two tables of this subsection. 

Annihilation neutralinos can also scatter off nucleons. For these we again take {E^) = 2Tc 
and ro = ^Rc, thereby obtaining the free-streaming conditions: 



M^(MeV) 





10 


30 


50 


100 


150 


200 


Mg(GeV) > 


154 


153 


147 


132 


152 


186 


214 



Thus, if Mg is sufficiently small the arguments of Sec. El cannot be used to set bounds on 
Me, since, below the values of Mg listed in this table, annihilation neutralinos will propagate 
diffusively out of the core because of strong neutron-neutralino scattering. 



6.2 Gravitational Trapping 

Interactions with matter are not the only way that neutralinos can be trapped in the supernova. 
Neutralinos produced in the outer regions whose kinetic energy obeys 

Eu. < (86) 

will be trapped by gravitational attraction. G is Newton's constant and R is the radius at which 
the neutralino is produced. Mji is the enclosed mass of the supernova at radius R. To get an 
estimate of the effect of gravitational trapping we take R ^ Rc and Mr ~ Msn- In terms of 
neutralino momenta, this means that neutralinos with |pj^| small enough that: 



will not escape from the supernova. Here 7 = l/Vl — v"^- Since GMgN / Rc ~ 0.15, this 
demonstrates that relativistic neutralinos (almost) always escape. 

In the non-relativistic case, it is easy to rewrite Eq. ()87|) as a condition on the neutralino 
velocity: 



< f , (88) 
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with Rs = 2GMsN ~ 4.1km the Schwarzschild radius of the core. Assuming that the distribu- 
tion of x velocities is Maxwehian ((i?kin) = 3Tc/2), we find for the non-relativistic neutrahnos: 

Combining this with Eq. ()88() we obtain a lower bound for the mass of gravitationally-bound 
neutralinos: 

> ( ^ ) 3Tc^ 285 MeV . (90) 



Rs 

This is outside the range of neutralino masses where we can set bounds {M-^ < 200 MeV) and 
thus our previous analysis holds. 



7 Conclusions 

In general in a supernova, neutralinos will be produced by both e~^e~ annihilation and NN 
binostrahlung. If the neutralinos free-stream out of the supernova, they can lead to excessive 
cooling which would alter the observed neurtino signal from SN1987A. In order to exclude this, 
we set bounds on the relevant supersymmetric masses: {M^, Mg) and (M^, Mq), respectively. 
The bounds are significantly stricter for e^e~ annihilation, since the production cross section is 
larger in that case. 

For e'^e~ annihilation we can, for a given neutralino mass, exclude the following values of 
selectron mass, Mg: 



M^(MeV) 





10 


30 


50 


100 


150 


200 


Me(GeV) < 


1275 


1275 


1260 


1188 


930 


700 


450 


Me(GeV) > 


322 


322 


321 


315 


341 


363 


377 



If the selectron mass is below the lower bound, the neutralinos produced by annihilation do 
not free stream, and so the arguments of this paper do not apply. Meanwhile, if Mg is above 
the upper bound then an insuffcient number of neutralinos are produced to affect the neutrino 
signal from SN1987A. 

For the case of the squarks we obtain a similar table 



M^(MeV) 





10 


30 


Mq{GeY) < 
Mg{GeV) > 


360 
306 


351 
305 


299 
298 



Only a small interval of Mq values is excluded, even for = 0, because there is only a narrow 
region where the coupling of the neutralinos to neutrons is large enough to compete with neutrino 
production, but small enough to avoid diffusive neutralino propagation. 

Furthermore, we expect that many-body effects will modify this squark-mass bound. Recent 
calculations of these effects in the case of neutrinostrahlung 29^ suggest that the dominant 
many-body correction to the binostrahlung rate at supernova temperatures will be the Landau- 
Pomeranchuk-Migdal (LPM) effect. The LPM effect will tend to suppress the production of 
axial radiation, thereby serving to still further weaken our already weak squark-mass bound [TTl 

isnmaiis]. 
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In Fig.^^we present the combined bounds for = 0. Supernova-cooling arguments exclude 
a large range of possible selectron masses if the neutralino is massless, but Mg > 1300 GeV is still 
allowed. In the region Mg < 300 GeV annihilation neutralinos diffuse due to electron scattering. 
They also diffuse, but this time due to interactions with neutrons, if Mq < 150 GeV. We cannot 
set bounds in these regions of parameter space. Their proper treatment would require a complete 
proto-neutro-star simulation. 
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Figure 11: Mg — Mq parameter space, with the regions excluded by the arguments of this paper 
for the case My. = shaded. The regions already explored by LEP searches are shown in grey. 
The white regions indicate areas of parameter space in which we are unable to set bounds. When 
this is because of diffusive propagation of the neutralino we have indicated the mechanism that 
leads to diffusion. Values of (Mg, Mq) in the upper-right-hand corner are permitted because there 
the neutralino production rate is too low to affect the SN1987A neutrino signal significantly. 

The bounds shown in Fig. ^Jdo not change significantly with M^. until masses of 50 MeV are 
reached. Above this value there is no squark mass bound whatsoever from SN1987A. Selectron 
mass bounds are also significantly reduced. For neutralinos with mass M^ > 200 MeV, we can 
essentially place no limits. Once this range of values is reached the LEP and Tevatron bounds 
on the sfermion masses, together with the Boltzmann suppression of heavy-particle production, 
guarantee that the supernova neutralino production rate is just too small to modify the neutrino 
signal markedly. 

In summary, we showed that a light neutralino, < 0(1 GeV) is allowed by labora- 
tory experiments and consistent with the MSSM if it is pure bino. We have investigated the 
bounds obtained from the SN1987A neutrino signal on a quasi-stable neutralino with mass 
M^ < 200 MeV and studied the production of such neutralinos in both electron-positron an- 
nihilation and nucleon-nucleon collisions. In the former case, we compared the Raffelt crite- 
rion with an estimate of the integrated emitted energy calculated using the radial temperature 
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and degeneracy profiles obtained in supernova simulations. We found good agreement between 
tlie two approaches. We then employed the Raffelt criterion to bound neutralino production 
in nucleon-nucleon collisions inside the supernova. Overall we find that for selectron masses 
300 GeV < Mg ^ 900 GeV we can exclude neutralino masses below 100 MeV. On the other 
hand, for selectron masses above 1200 GeV there is no lower bound on the lightest neutralino 
mass. There is also a near absence of any bound on the squark masses: only a narrow range of 
Mq can be excluded, and this only for < 20 MeV. 
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A Appendix 

In this appendix we give the explicit result for the angle-averaged hadronic tensor. Tin, in terms 
of phase shifts. The quantity needed for the leading order diagrams is the nucleon-nucleon 
scattering T-matrix T/vat in plane wave decomposition. It can be decomposed into partial 
waves in the following way. (We only give the expressions for equal final and initial spin. This 
is all that is required for the reaction nn — > nnxx-) 



{SM's,p'\Tnn\SMs,p) = Y,YL'M',{p')YLMAPr'^(j^,S'M's, S Ms, L', AL,Mj) , (91) 



Tip,S'M's, SMs, L', AL,Mj) = ^,{3' M's, L' Mi\J Mj){S Ms, LMl\JMj) T^^^'^^'^ip) . 



where 5^ denotes the phase shift in the partial wave denoted by L These quantities are available 
online Note, only elastic scattering is kinematically allowed under the conditions of the 

study presented. 

It is straight forward yet lengthy to relate T to the angle-averaged hadronic tensor TYjj as 
defined in Eq. (f7T|) : 



where 





(93) 



Y,'^{{2-MsM's)\T[p,lM's, I Ms, L', AL,Mj)\^ 
-2Re [T{p, l{M's-l),lMs, L', AL, Mj - 1)* 

xT{p, 1 M's, 1 {Ms + 1), L', AL, Mj)] } . 



(94) 
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